Quench dynamics near a quantum critical point: application to the sine-Gordon model 
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We discuss the quench dynamics near a quantum critical point focusing on the sine-Gordon model 
as a primary example. We suggest a unified approach to sudden and slow quenches, where the tuning 
parameter \{t) changes in time as A(t) ~ vf^ , based on the adiabatic expansion of the excitation 
probability in powers of v. We show that the universal scaling of the excitation probability can be 
understood through the singularity of the generalized adiabatic susceptibility X2r+2(A), which for 
sudden quenches (r = 0) reduces to the fidelity susceptibility. In turn this class of susceptibilities 
is expressed through the moments of the connected correlation function of the quench operator. 
We analyze the excitations created after a sudden quench of the cosine potential using a combined 
approach of form-factors expansion and conformal perturbation theory for the low-energy and high- 
energy sector respectively. We find the general scaling laws for the probability of exciting the system, 
the density of excited quasiparticles, the entropy and the heat generated after the quench. In the two 
limits where the sine-Gordon model maps to hard core bosons and free massive fermions we provide 
the exact solutions for the quench dynamics and discuss the finite temperature generalizations. 



; I. INTRODUCTION 

The Sine-Gordon (SG) model is one of the few exam- 
ples of exactly solvable models both in the classical and in 
the quantum case. In the classical limit, an arbitrary ini- 
tial perturbation splits into solitons and breathers, which 

' are infinitely long lived excitations (see e.g. Ref. 

In the quantum regime these solitons and breathers are 
quantized and form a discrete spectrum (for each mo- 
mentum state) 0, 0] . This model found numerous ap- 
plications in various fields of physics including statistical 
physics, condensed-matter physics, atomic physics and 
high-energy physics. Its integrability allows to solve var- 
ious equilibrium and non-equilibrium situations and gain 
important insights into phenomena ranging from classical 

] and quantum phase transitions [^-Q to quench dynamics 
of coupled superfluids 0, [l] . 

\ For the purposes of this work we will keep in mind the 
applications of this model to one-dimensional quantum 
systems. Thus a natural low-energy theory describing a 
one-dimensional interacting gas of bosons or fermions is 
the Luttinger liquid model, which is nothing but a scalar 
free 1-1-1 dimensional bosonic theory [5(1 . The Luttinger 
liquid is often unstable to various perturbations opening 

\ a gap in the system coming either from the external po- 
tential (periodic Q or disordered [10]) or the tunneling 
coupling between multiple Luttinger liquids [ll| . In these 
situations the most relevant perturbations appear in the 
form of the lowest harmonic of the cosine potential al- 
lowed by the symmetry [6l.lT^ leading to the sine- Gordon 
Hamiltonian: 

H=^Jdx [{dMf + (n(a;))2 - 4Acos(/3</)(a;))] . 

(1) 

In cold- atom systems this Hamiltonian can be imple- 
mented through optical lattices [9,, JJj, it also emerges 
naturally in the context of interacting quantum wires 
@ : [3 ' mesoscopic superconducting junctions [l^ , spin 



chains 0] and many others. 

The recent experimental developments in cold-atoms 
[igI [17| prompted a rapidly growing theoretical interest 
to study the dynamics of quantum systems. A particu- 
lar attention was paid to such issues as sudden quench 
dynamics [l8l - [23 |. thermalization in integrable and non- 
integrable systems (25l - l28j and slow dynamics in isolated 
systems near quantum critical points (29l-l3TI|. The latter 
works illustrated the universal scaling of the density of 
quasiparticles generated due to nonadiabatic transitions 
during the crossing of quantum phase transitions. The 
universal properties of quantum critical points in equi- 
librium extend to the dynamics and allow one to make 
universal predictions on the defects and energy genera- 
tion similar to the Kibble-Zurek relation s fs^lssl for slow 
dynamics near quantum critical points |34l - l44l | . 

In the preceding work [isj we showed that a similar 
universality also persists for a sudden quench dynamics 
near quantum critical points if the quench amplitude is 
sufficiently small. There we analyzed the general scal- 
ing behavior of different quantities following a sudden 
quench of the amplitude A/ starting from the critical 
point Xi — 0. We presented some general arguments 
for the scaling of various quantities such as the proba- 
bility of exciting the system Pex, the density of excited 
quasiparticles ricx, the (diagonal) entropy density Sd [i^ . 
and the heat density Q [47[, with the quench amplitude 
and the system size. We showed that the correspond- 
ing scaling laws are universal and can be described by 
two critical exponents characterizing the quantum criti- 
cal point: the dynamical exponent z and the correlation 
length exponent i'. We also argued that the scaling be- 
havior for sudden quenches is smoothly connected to the 
scaling behavior of similar quantities for slow quenches, 
where the tuni ng par ameter is turned on as an arbitrary 
power of time [4^, ll^: A(i) ~ vt^, where u is a small pa- 
rameter. In particular, in the limit r — this parameter 
becomes the quench amplitude: u = A/. 

The goal of this paper is to verify how the general ar- 
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guments of Ref. [45| work in practice. We will use the 
sine-Gordon model as a primary example. Performing 
explicit calculations based on the adiabatic perturbation 
theory, form-factor approach, and conformal perturba- 
tion theory we will illustrate how these scalings are re- 
alized for the case of quenching the cosine potential and 
highlight various subtleties which can emerge due to pos- 
sible ultraviolet or infrared divergencies. In the two lim- 
its where the SG model reduces to free massive bosons 
and fermions we will extend the analysis to higher dimen- 
sions and again verify the general predictions of Ref. (45| . 
In these two limits we will be able to generalize the re- 
sults to finite temperatures and see how the quasiparticle 
statistics affects the scaling of various quantities. 

The paper is organized as follows: in Sec. HI] we re- 
view the general scaling results presented in Ref. 45] in- 
troducing a generalized framework that allows us to de- 
scribe sudden quenches and slow quenches on the same 
footing. In Sec. IIIII we describe the basic properties of 
the SG model. Then in Sees. IIVI and |V] we analyze the 
quench dynamics in the SG model using the adiabatic 
perturbation theory. In Sec. I VII we present the exact 
analysis of quenches in the two noninteracting limits of 
the SG model corresponding to the free massive bosons 
and free massive fermions. In Sec. IVIDI we furthermore 
extend this analysis to finite temperature quenches. Fi- 
nally in Sec. I VIII we discuss some interesting relations 
between quench dynamics, fidelity susceptibility and ge- 
ometric tensors. 



II. GENERALIZED FRAMEWORK FOR THE 
ANALYSIS OF QUENCH DYNAMICS: 
HIERARCHY OF SUSCEPTIBILITIES 

In this section we will illustrate a generalized approach 
to study the scaling laws of the quantities of interest after 
performing a quench. As in Ref. (45j we consider a d- 
dimensional system described by a Hamiltonian H{\) — 



Hq + XV, where Hg is the Hamiltonian corresponding to 
a quantum critical point (QCP) and ^ is a relevant (or 
marginal) perturbation. In particular, we consider the 
processes where the coupling changes near the critical 
point as a power law in time: 

x{t) = v^^eit), (2) 

r! 

where u is a small parameter that we define as v = 
\t=o and 8 is the step function, so that the cou- 
pling starts changing with time as near the critical 
point. On the same footing we can consider the opposite 
processes, where quenches stop at the critical point as: 

m^v^^i^eitf-t). (3) 

The analysis also extends to cyclic processes where the 
dynamical evolution starts and ends at the quantum crit- 
ical point. E.g. X{t) X vf{tf -ty/irlfj). In all these 
cases V plays the role of the adiabatic parameter such 
that u — ^ corresponds to the adiabatic limit. In par- 
ticular, if the system is initially prepared in the ground 
state then in the limit u — > it will remain in the ground 
state. The physical meaning of v is quite transparent. 
For sudden quenches r = 0, it corresponds to the quench 
amplitude f = A/; for linear quenches r = 1, it describes 
the quench rate v = X; for quadratic quenches r = 2, 
it describes the acceleration near the quantum critical 
point: V — X and so. The smallness of v suggests the 
possibility of using it as an expansion parameter. Rely- 
ing on the adiabatic perturbation theory one can indeed 
write an expansion in the powers of v. The details of 
this approach can be found in Refs. [13, IMI- We only 
quote the final expression for the transition amplitude to 
the excited state |n) as a result of this process (see e.g. 
Eq. (19) in Ref. [50] or Eq. (18) in Ref. [HI): 



an{tf) 



iX 



{n\dt\0) 



{n\m 



En [t) - Eo (t) En (t) - Eo it) dt En {t)-Eo it) 
n\dx\0) V {n\d^\0) .-o 1 



A- 



EniX)^EoiX) iEniX) - EoiX)) 



X' 



{n\dx\0) 



En{X)-EoiX) dX EniX) - EoiX) 



_»(e„(A)-eo(A)) 



(4) 



where &nit) = EniT)dT is the dynamical phase (in 
general one should also add the Berry phase term) and 
En is the eigenenergy of the state \n) . We can anticipate 
that the critical point will dominate the dynamics so that 
only one of the limits A = Ai (if Ai = 0) or A = A/ (if 
A/ = 0) should dominate the transition amplitude. From 



this equation it is clear that the lowest non-vanishing 
time derivative of A(i) at the critical point dominates 
the asymptotics of the transition probability: 

I ,2^ 2 \{n\dxm' \{n\vm' 

(EniO) ~ EoiO))^^ (i?„(0)-i?o(0))2'-+2' 

(5) 
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where all the matrix elements are evaluated at the quan- 
tum critical point. Therefore we find that the probability 
of excitations, defined with respect to the instantaneous 
ground state, scales as: 



Pcx = ^|a„P«iVx2,.+2(0), 

n#0 



where 



\{n\V\OW 



(6) 



(7) 



is the generalized adiabatic susceptibility of order m. For 
m — 2 (r = 0) X2 = X/i where x/ is the fidelity suscep- 
tibility: 



Xf 



1 
I? 



mV\n)\' 



(8) 

Thus Xf describes the system's response to a sudden 
quench (see Refs.ji^ and [l^). For m = 4 {r — 1) the 
linear adiabatic susceptibility X4 describes the response 
to a linear quench and so on. Let us note that the scal- 
ing dimension of the susceptibilities Xm(A) is related to 
the scaling dimension of the fidelity susceptibility via: 
dim[xm(A)] = dim[x/(A)] — z(m — 2) — d—2/i^ — z{m — 2). 
If dim[xm(A)] becomes negative {di' < 2+zv{ra—2)) then 
the susceptibility Xm(A) diverges at the quantum critical 
point. Since dim[A] — 1/v the asymptotical behavior of 
Xm(A) at small A is: 



Xm(A) 



(9) 



In finite size systems this divergence is cutoff by the sys- 
tem size L when the correlation length ^(A) ~ l/IA]" 
becomes comparable to L. Thus exactly at the quan- 
tum critical point the generalized adiabatic susceptibility 
scales as 



Xm(0) 



(10) 



For ra — 2 this scaling agrees with that of the fidelity 
susceptibility jssl . Is^ ] . If the scaling dimension of Xm is 
positive then this generically indicates that this suscep- 
tibility has a cusp singularity, meaning that the asymp- 
totics (HI) and (ITOl) become subleading. 

Combining Eqs. © and pU)) we find that if the condi- 
tion dv < 2{1+ zvr) is fulfilled, then Pex has the following 
asymptotic behavior in the adiabatic limit: 



2j2/p+2zr 



(11) 



Since the probability Pox must always be smaller than 
one, this scaling must be valid only for quenches of small 
amplitude |u| < 1/^1/"+^^ The physical meaning of 
this requirement is that the correlation length of the sys- 
tem associated with this velocity ^(u) ^ l/|u| must 



be big compared to the system size ^ ^ L. Physically 
the length scale ^{v) characterizes the crossover in the 
response of the system from the sudden to the adiabatic 
regimes similar to the Kibble- Zurek mechanism [sH, fs^. 
I.e. for excitations with characteristic size (e.g. wave- 
length) shorter than ^(u) the quenching process looks 
adiabatic, while for excitations with size larger than £,{v) 
it looks sudden. For linear quenches ^(u) ^ in- 
deed reproduces the characteristic Kibble-Zurek length 
scale [32l.[33j. 

Let us point that the crossover to the regime with 
^ex(f) ~ 1 occurs precisely when the characteristic 
length scale ^(u) becomes of the order of the system size. 
This can be explained using a simple symmetry argu- 
ment, which is easier to understand for sudden quenches 
r — 0. In this case the quench corresponds to projecting 
the ground state in the critical state on the new basis. 
The quantity 1 — Pcx defines the probability that the sys- 
tem remains in the ground state after the quench. But if 
^(u) < L the new ground state has a well defined sym- 
metry absent at the critical point. Thus the probability 
to remain in the ground state should be close to zero. 
Similar considerations apply to slow quenches, r > 0. 
Thus for quenches of higher rate \v\ 3> l/L^^'^'*'^'' (but 
still slow compared to the other microscopic scales of the 
system) the probability of exciting the system is close 
to one and is no longer a good quantity to character- 
ize the system response (like the overlap of two ground 
states becomes very close to zero and thus not very infor- 
mative). We argued that in this case the physically 
relevant and meaningful quantity to look at is the density 
of generated excitations riox (quasiparticles if we are deal- 
ing with nearly integrable models). Indeed the physical 
reason why the global probability of excitation rapidly 
approaches unity is that exciting even a single quasipar- 
ticle makes the new state orthogonal to the initial state. 
Yet physically we anticipate that exciting a single quasi- 
particle in the system should not change drastically the 
properties of the system. In particular, it should not af- 
fect the probability to excite the next quasiparticle. Since 
we are typically dealing with few body operators by per- 
forming a quench we couple only to the many-body states 
characterized by few quasiparticles. Thus it was argued 
that the density of generated quasiparticles should scale 
as: 

for \v\ < l/Xi/'^+z'- and 

ricx - (13) 

in the opposite limit. Note that the second scaling has 
a much wider domain of applicability than the first one 
since there is no requirement on vanishing quench ampli- 
tude with the system size. In particular, for r — 1 the 
result ((T5)) agrees with the Kibble-Zurek scaling predicted 
earlier in Refs. [2^|33|. Also for nonlinear quenches r > 1 
the scaling agrees with the results of Refs. [1^ l49| . 
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noting that v = (1/r)'" in the notations of Ref. [48| . 
Let us also point out that, as we wih ihustrate here for 
the specific case of the sine-Gordon model (Sec. |Vl, for 
dv > 2 the nonanalytic asymptotic with |u| or L in the 
scahng for and Pcx survives but becomes subleading. 
The leading asymptotic is then analytic and it comes 
from the excitations of non-universal high energy states. 

The density of excitations is a very simple and intu- 
itive object, however it is not always well defined. Un- 
less we deal with long lived quasiparticles or other e.g. 
topological excitations, in nonintegrable systems riox is 
not conserved in time after the quench because of var- 
ious relaxation processes. For the same reason ncx is 
not a readily observable quantity. Hence one needs to 
characterize the response of the system by other means. 
The two natural quantities, which can be defined for any 
Hamiltonian system, are the (diagonal) entropy [46| and 
the heat (or the non-adiabatic part of the energy, which 
is also the excess energy above the new ground state of 
the quenched Hamiltonian [13]). Because both quantities 
are extensive, it is convenient to deal with their densities: 



(14) 



Sd - --^^|a„plog 



n n^O 

Here the index n enumerates the many-body eigenstates 
of the Hamiltonian in the final state and |q;„P, as written 
above in Eqs. (|H) and (O), is the probability of the tran- 
sition to the n-state due to the quench. The scaling of the 
entropy is very similar (up to possible logarithmic correc- 
tions) to the one of Pex and nex- The advantage of the 
entropy over Pex is that it is well behaved and meaning- 
ful for large amplitude quenches 3> l/L^/""'"^'', where 
-Pox 1. The scaling of the heat is different from that 
of Ticx because of the term En — Eq in Eq. (fT5|) . If the 
quench ends at the critical point then one finds that the 



scaling of the heat is described by X2r+i- 
so that 

for \v\ < l/Li/''+^'' and 



Q 



(16) 



(17) 



for \v\ > l/Li/''+^''. 

Let us note that the susceptibilities Xm are directly 
related to the connected parts of the correlation func- 
tions Ga(t) = {V{X, t)V{X, 0)) - {V{, A, 0))2 by a simple 
generalization of the expressions obtained in Ref. [53] : 



1 



U\m~ 1)! io 



^Gx{T)dT. (18) 



For local perturbations V — Vi^) it is obvious that 
if G{t) decays sufficiently fast at large r then Xm is a 
non-negative number well defined in the thermodynamic 



limit L — > oo. For the sine-Gordon model we wrote the 
exphcit expressions of xi and X2 in Eqs. ([5^ and ([55]) . 
Using the Lehmann's representation of the correlation 
function we can define the generating functions for the 
susceptibilities Xm(A): 



^ ^ \{n\Vm^ 
rd 



^^En{X)-EoiX) + a 



such that 



/(a) 



i-arxM- 



(19) 



(20) 



This characteristic function clearly contains the informa- 
tion about the scalings laws of the quantities discussed 
in this work for all possible power law quenches. 

Let us say a few words about the scaling in the regime 
where the relevant adiabatic susceptibility has only a 
cusp singularity. In this case the leading asymptotics 
of Xm(A) at A — >■ is given by a nonuniversal constant, 
which in turn defines the leading response of the system 
to the quench. Thus one can expect the scaling for the 
quantities n^x, S^, Q to be analytic in v. Instead of 
Eq. (fl^ we anticipate: 



"ox ^ t^^X2r+2(0). 



(21) 



Similarly for the heat (in the case of a quench ending 
at the QCP) we expect Q ~ v'^X2r+i{0)- As we will 
show below the non-analytic behavior of these quantities 
with V survives, but becomes subleading and thus hardly 
identifiable either numerically or experimentally. 

It is interesting to note that these scaling results ex- 
tend to the situation of gapless systems even in the ab- 
sence of the critical point. Thus if we are dealing with 
gapless bosonic theory then one simply needs to take the 
limit — >■ CX) in the expressions above. It is clear from 
Eqs. ([9|) and (fTO|) that for m > 2 the susceptibility Xm 
can still diverge in low enough dimensions if d < z{m—2). 
This divergence is actually the origin of the strong non- 
adiabaticity in low-dimensional gapless systems [55] (see 
also Ref. [5a]). Thus the existence of the quantum criti- 
cality is not really needed to have a non-analytic scaling 
of various observables with v. Following Ref. [5^ let us 
make the following observation. In the limit z/ — > oo in 
low dimensions the scaling of the three quantities riex, Sd 
and Q is only sensitive to the quench time but not to the 
actual quenching protocol. For instance for d < z{2r— 1) 

the scaling of the heat is Q ~ ju]^'^^ ~ l/r^~, where 
T ~ (A//u)^/'" is the quench time required to change A 
from the initial value A^ = to some fixed final value A / . 
Thus the scaling of Q with r is not-sensitive to the power 
r with which the tuning parameter is quenched. On the 
other hand for d > z(2r — 1) the situation is opposite. 
The scaling of the heat Q ^ ^ I/t^"^ is very sensi- 
tive to the power with which the coupling A is turned on 
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(off). This suggests that the optimum quench protocol, 
which gives the smahest excitation level in the system 
for a given quench time r, corresponds to the crossover 
power r = ^(f + l)- Higher powers of r are to be avoided 
because of the growth of the prefactor with r j49|] . This 
argument extends the predictions of Ref. [i^ to generic 
gapless systems. 



III. THE SINE-GORDON MODEL 



The SG-model is described by the Hamiltonian: 



H=\l dx [{dA{x)f 



(n(a;))2 -4Acos(/3(/.(a;))] , 



(22) 

where 0(a;) and n(a;) are conjugated fields. Since this 
model is exactly solvable a great amount of facts is known 
about it. It is conventional to introduce the parameter ^ 

M 



(23) 



The spectrum of the SG Hamiltonian ([2^ depends on 
the value of (5 (or i^). We will be primarily interested 
in the regime of < /^^ < Stt since in this range for a 
finite A the system is in a gapped phase, while for A = 
it corresponds to the gapless Luttinger liquid. Therefore 
the system in this regime undergoes a quantum phase 
transition at A = 0. Furthermore for Att < 0^ < Stt 
(1 < f < oo), so called repulsive regime, the spectrum of 
the SG Hamiltonian consists of solitons and antisolitons. 
At fixed small A and /S^ = Stt the systems undergoes 
a Kosterlitz-Thouless transition to the Luttinger liquid 
regime, where the cosine term is irrelevant. The point 
0^ = 47r (^ = 1), known as Toulouse point, maps to 
the free massive fermionic theory; at this point solitons 
and antisolitons correspond to particles and holes. For 
< /J^ < 47r (0 < ^ < 1), so called attractive regime, the 
spectrum in addition to the solitons and antisolitons con- 
tains their bound states called breathers. The number of 
different types of breathers depends on the interaction 
parameter ^ and is equal to the integer part of 1/^. We 
denote a breather by i?„ with n = 1, 2, . . . , [1/^]. In 
the small /3 (or ^) limit the SG model is well described 
by the gaussian approximation, where one expands the 
cos(/3(/)) term in the Hamiltonian (22) to the quadratic 
order in </>. In this limit there is only one massive excita- 
tion, corresponding to the lowest breather Bi. Solitons 
(kinks), antisolitons (antikinks), and breathers are mas- 
sive particle-like excitations. The soliton and antisoliton 
mass in terms of the parameters of the Hamiltonian ((22)) 
was computed in Ref. 0]: 



^r(^ A\ 2r(| 



(24) 



Breathers mass is related to the soliton mass via: 
Mb =21/3 sin' 



(25) 



Note that for weak interactions (small ^) the lowest 
breather masses are approximately equidistant, Mb^ oc 
n, which suggests that these masses correspond to eigen- 
energies of a harmonic theory. There is a direct analogy 
between the breathers in the sine-Gordon model and the 
energy levels of a simple Josephson junction Q. In the 
Josephson junction (which is described by the Hamilto- 
nian (P^l without the first gradient term) the energy lev- 
els also become approximately equidistant if the interac- 
tion (charging) energy is small. Let us note that because 
the SG theory is explicitly relativistic invariant the dy- 
namical exponent z is equal to 1. The soliton mass 
identifies the characteristic energy scale of the system, 
therefore we expect for it a universal behavior around 
the critical point: Ms ^ lA]^'', that defines the critical 
exponent v for this theory (from Eq. ([M]) ) 

^-'-^ (26) 

or V = 47r/(87r — /3^), which can be continuously varied 
between 1/2 and 00. 

We will focus on sudden quenches of the parameter A, 
A(t) = A/0(t) and we will illustrate how the general scal- 
ings laws for the density of excited quasiparticles (Eqs. 
(fT2l) and (fT3)) ) can be obtained for this model. As we dis- 
cussed in the previous section and in Ref. these scal- 
ing laws are closely related to the scaling of the fidelity 
susceptibility Xf- This susceptibility can be expressed as 
an integral of the correlation functions. For the SG model 
the latter can be accessed from two limits: the infrared 
limit (IR) describing long distance properties, and the 
ultraviolet limit (UV) corresponding to short distances. 
The approach based on the form factors is well suited to 
analyze the former IR limit. To deal with the UV limit 
one can rely on the conformal perturbation theory, which 
is based on the fact that at short distances the SG model 
(as well as many other IR-massive models) is free and 
conformal (it effectively reduces to the Luttinger Liquid 
theory). These two approaches are complementary and 
can be used depending on the physics of interest. As 
we will show below for dv < 2 (or v < 2 since d = 1) 
the former form factors approach reproduce the correct 
scalings of Pox and riox, while ioT v > 2 one should use 
the conformal perturbation theory. Similar story is valid 
for heat but the IR and UV dominated domains are de- 
fined according to whether (d + z)v = 2v is bigger or 
smaller than 2. In general, it is rather difficult to relate 
the results from the two formalisms: in order to get the 
UV physics from the form factors expansions one has to 
sum up the large (if not infinite) number of contributions 
(see Ref. Q and the following Sec. |lV]and|V|. The con- 
verse is also true. We will analyze both the UV and IR 
contributions to the excitation probability and show that 
their combination indeed reproduces the correct scaling 
asymptotics. 
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IV. ANALYSIS OF THE LOW ENERGY 
EXCITATIONS BASED ON THE FORM FACTOR 
EXPANSION 

Let us start from the perturbative expression for the 
probabiUty of exciting the system based on the adiabatic 
perturbation theory. Within the latter the transition am- 
phtude to the state \n) is found as [50| 



d\{n\dx\0). 



(27) 



If the matrix element (n|(9;j|0) is not singular at the QCP 
then this expression reduces to the result of the first order 
of the conventional perturbation theory. The expression 
(l?7l) has an advantage that it symmetrically treats both 
initial and the final values of the coupling A and gives 
convergent expressions even when the conventional per- 
turbation theory breaks down. From Eq. ([27l) we find 



E 



d\id\2{0\dM{n\dxM- (28) 



In general one has to sum over all intermediate states \n). 
However, as we argued in [45], \i dv < 2 (i.e. ^ < 3 in 
the SG model) then Pox is dominated by the low energy 
excitations. Since solitons and breathers are massive, i.e. 
gapped for |A| > 0, we can expect that the dominant con- 
tribution to Eq. (pS)) comes from single pairs of the lowest 
energy excitations (solitons and antisolitons for ^ > 1 and 
Bi breathers for ^ < 1) with opposite momenta. There- 
fore the probability of exciting the system and the num- 
ber of excited quasiparticles will become identical (up to 
a factor of two reflecting that one excited state creates 
two solitons or breathers). In this case we can identify 
the sum over the states n as a sum over momenta k and 



use: 



afc(A/) 





A/ 



2 / d\ 



d\'{k\dy\Q) 

, (O|cos(/30)|fc) 



(29) 



where \k) is a short-hand notation corresponding to the 
soliton-antisoliton (breather) pair with momenta k and 
— k and -Bfe(A) is the energy of such a pair. The expression 
for the density of excited quasiparticles riex is therefore: 



1 --^ /■ dk P'^f 

(A/) « ;^El"'=(^/)l'= J ^ I dX'{0\dy\k) 



(O|cos(/30)|fc) 



(30) 



Note that a factor of 2 should appear here, due to the 
fact the excited states correspond to pair of excitations. 



but it is exactly canceled since we are counting each pair 
twice integrating over all positive and negative momenta. 
Likewise in the same order of approximation one can ob- 
tain the expressions for the heat density Q: 



dk 



^' (O|cos(^</>)|fc) 
^ Ek{X') - Eo{X') 



and the entropy density: 

„ 1 rdk. 



2 / 27r 



(31) 
(32) 



There is an additional factor of 1/2 in the expression 
for the entropy coming from the fact that particles are 
excited only in pairs and the contribution of each pair to 
the entropy (unlike ricx) is not doubled. 

As we argued earlier (see also Ref. [i^) the advan- 
tage of using intensive quantities like riex, S'd, Q over Pgx 
is that the regime of validity of the perturbation theory 
for them is much bigger. Unlike for Pcx the quench am- 
plitude is not required to vanish as some power of the 
system size. In Sec. I VII we will explicitly illustrate this 
point for two exactly solvable limits of the SG model: 

^ = 1 and C < 1- 

The adiabatic perturbation theory allows us to reduce 
the computation of the dynamical response to static cor- 
relation functions. The same is true if we are dealing with 
a slow processes where the coupling changes gradually in 
time. In this case (see Sec. [ll]) there is an additional 
dynamical phase entering the transition amplitude. For 
example for linear quenches X{t)= vt instead of Eq. ((29)) 
we should use (see also Refs. pfsl, Isol, [57|): 



ak{v) = / dX 



^, (O|cos(/?0)|fc) 
Ek{X') - EoiX') 



X exp 



A' 



dX" {Ek{X")-Eo{X"}) 



(33) 



in all expressions for riex, Q and Sd- 

The matrix elements appearing in Eqs. ([^ and ([55)) 
are related to the form factors of the operator cos(/30). In 
general the form factors represent the matrix element of a 
particular operator between the vacuum and the asymp- 
totic states (eigenstates of the Bethe ansatz) created by 
the Zamolodchikov-Faddeev operators corresponding to 
solitons, antisolitons and breathers. Because the SG the- 
ory is Lorentz-invariant it is convenient to use the rapid- 
ity variable 0, — oo < 9 < oo, which parametrizes the 
energy and momentum of the soliton (breather): 



E = Ms,B„ cosh e, 



k = AL 



sinhfl. 



(34) 



The eigenstates of the SG model can then be labeled by 
\0n ■ • ■fi'i)a„...ai, whcrc fl^ = {s, s, n = 1 . . . [1/C]} corre- 
spond to solitons, antisolitons, and breathers. The cor- 
relation functions of an arbitrary operator 0{x,t) are 
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written as an infinite series expansion in terms of all the 
asymptotic states j58| : 

(o(.,<)o.(o,o))=f;i:/np^ 

e'E-.ifc^.-iJ,* i^Op^Q^O)!^^^ . . .0i)a„...aj' .(35) 



The matrix elements entering Eqs. ((29)) and ([SSj) are 
then explicitly related to the form factors of the oper- 
ator cos(/30) in the two-particle asymptotic states: 



|cos(/30)|fc) = 



1 



cos(/3./))|0,-0)),,,a,. (36) 



Here the additional factor of XjEiff) comes from the rela- 
tivistic normalization of the eigenstates parametrized by 
6», {Q\Q') = S{0 - 9')/{2n). While in Eq. §UI we used a 
different normalization {k\k') = 5{k — fc')/(27r). Taking 
into account the relations (|34|) we indeed find that the 
normalization factor l/\dgk{9)\ — 1/E{9). 

Let us discuss the two-soliton contribution to ricx first, 
which is expected to be dominant for 1 < £ < 3. The 
corresponding form factors in the SG model are known 
in the literature [H^, HO] : 



cosh(6l/2)exp[o-^] 
sinh(^) ' 



(37) 



where we kept only the relevant factors (those which de- 
pend on either 9 or A/). We point out that the sum- 
mation over a comes from considering both the possi- 
bilities for a pair of excitations (soliton-antisoliton and 
antisoliton-soliton) . Also we write both signs in the ex- 
ponent of the operator e^^^'^ since these form factors are 
invariant with respect to the change of sign. Furthermore 
it is known that: 



5^ 



(38) 



This result together with Eq. implies that 

(O|cos(/30)|fc) - \Ms\^^P{9) ~ \\f\^P{9), where 
P{9) is some function which only depends on 9. Next 
we use this scaling of the matrix element in Eq. ([50)1 
for finding riox- If the integral over k (or equivalently 
over 9) converges at small k ^ Mg then we can perform 
the appropriate rescaling of the variables (A' — > A'/A^ 
and k — ?> k/Ms and send the upper limit of integration 
over fc to oo. Then we immediately obtain the scaling 
TT-ox ^ |Afs| ~ I'^/Tj which agrees with the general pre- 
diction (|13p in the case of a sudden quench. A simi- 
lar analysis of the matrix elements leads to the correct 
scaling for slow quenches characterized by an arbitrary 
exponent r (see Ref. [lH for linear quenches). 

To analyze the convergence of the integrals we need to 
find the asymptotical behavior of the matrix element in 
Eq. (15(11) at large 9. In order to do this, we look at the 



large 9 behavior of the ^-dependent piece in Eq. ([37]). In 
particular we find: 



G[29) exp 



±- 



1 /I 



"2 



1 



±oo, 



(39) 



which could be derived from the infinite-product-gamma- 
fimctions representation of G{9) 61] (or also, in a differ- 
ent way, see Eq. (61) in the book by F. Smirnov [59]). It 
is also straightforward to deduce the asymptotics of the 
remaining part of the form factor to finally get: 



^COS(^0) 

(29) 



exp 



1 

2e 



(40) 



To find the matrix element (0) cos(/?0)]fc)) we need to 
divide this asymptotical form by the soliton energy (see 
Eq. dMl) E = M,cosh(6i) w Msexp[6']/2 « fc at large 9 
so that 



(Olcos(^0)lfc) - M, 



ITT 



exp 



2 £ 



M. 



5-1 



E 



(41) 

This result correctly reproduces the scaling in the free 
fermionic limit of the theory (£ = 1), which can be ob- 
tained by elementary methods (see Sec. I VI C]) . Such scal- 
ing implies a logarithmic singularity at ^ — I for the 
heat density in a sudden quench Q ^ jA/J^logJA/] (see 
Eq. (pij) . which is equivalent to the log-divergence of the 
ground state energy. This divergence is expected because 
for £ = 1 we have v = {l+£,)/2 = 1, which corresponds to 
the crossover point between the scaling Q ^ ]A/1(''+^^'^, 
expected for v < 2/{d + z), and the quadratic scaling, 
expected for ly > 2/{d + z). For £ > 1 the integral over 
k in Eq. pip diverges at large fc, i.e. the main con- 
tribution to the heat density comes from high energies 
(UV limit), where the lowest form factors expansion is 
not justified and we need to rely on the UV limit of the 
matrix elements discussed below in Sec. |Vl As for the 
other three quantities: nox, ^cx, and Sd the scaling of 
the matrix element (PT|) together with Eqs. (15(11) and ([5^ 

suggests that riox ^ j'^/j'''' ~ \-^f\~^ holds for the whole 
soliton region 1 < £ < oo and there is no indication of 
crossover to the quadratic scaling at £ = 3 (correspond- 
ing to = 2). This result is however misleading because 
it comes from the fact that the lowest order form factor 
expansion does not correctly reproduce the ultra-violet 
asymptotics of the matrix elements. As we mentioned 
above to get the right UV limit one needs to perform an 
infinite re-summation of the form factors which implies 
that the sum in Eq. (1281) should be taken over multiple 
soliton states. This problem is well known in the equilib- 
rium analysis of the UV asymptotics of the static correla- 
tion functions ^62,] . In the next section we will show how 
the power two emerges from analyzing Eq. (|28p using the 
conformal perturbation theory. 

Likewise one can analyze the scaling of the matrix el- 
ements in the breather region (£ < 1). Since we have 
z/ = (1 -t- £)/2 < 1 we expect that the scaling of all the 
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quantities we are interested in Q , P^^, ric^, S d is domi- 
nated by low energies where the lowest form factor ex- 
pansion is justified. So for ^ < 1 we need to repeat the 
analysis above with the only difference that the states 
\k) in Eq. p9|) denote pairs of breather states with op- 
posite momenta. The form factors for breathers have a 
different scaling form than those for solitons. With the 
operator cos(/30) it is possible to have non-zero form fac- 
tors of a single breather B„ of the form (0| cos{(3(j))\Bn) 
if n is even. Those form factors do not depend on the 
rapidity at all (see [6^ for more details). The first 0- 
dependent contribution into Eq. (|28p comes from exciting 
two Bi breathers. These breathers are the only addition 
to the solitons in the spectrum for 1/2 < ^ < 1 and 
they are also the only surviving excitations in the non- 
interacting bosonic limit of the SG model ^ — 0. There- 
fore here we will concentrate only on the scaling anal- 
ysis of ricx and other observables coming from exciting 



pairs of Bi breathers with opposite momenta [821]. The 
matrix element describing this process is found through 
the corresponding form factors which were computed in 
Ref. [63| using bootstrap, the procedure of fusion of sev- 
eral breathers [6^] : 



(O|cos(/30)|fc) = 



where 



2 cos 



2 



and 



1 



E{e) 



(O|cos/30|Bi(0)Bi(-0)) 



1 



2Af, cosh 



/2sin 



2 



exp 



''^ dt t 



2tt sint 



(43) 



R{0) = exp 



dt sinh(t) sinh(t^) sinh(t(l + ^)) 
t sinh^(2i) 



exp 



dt sinh(t) sinh(iO sinh(t(l + 0) • , 2 / 

^ smh t { 1 — i — 

t sinh2(2<) 



r 



(44) 



Since the dependence of the matrix element (H^ on Mg 
(and hence on A/) is again found via the function Qp, 
we conclude that the scaling of n^x (as well as Q, Sd and 
Pcx) on \f is the same as in the repulsive regime, e.g. 
"■Gx ^ l-^/l'"^^^''^^- We only need to verify that the inte- 
grals over momenta k converge at large k. In order to 
do this we need to know the asymptotics of the function 
R{26) at large 9. The scaling at large rapidity for the 
function R{9) (defined in ^]) is R{0) 

^ ^Bi ' which does 

not depend on 9. This saturation, which is the conse- 
quence of bootstrap, is actually true for higher order form 
factors as well. For example (0| cos/30|i?2(6')-B2(— ^?)) ^ 
{Fb^Y ■ On the other hand, because of the bootstrap rela- 
tion, this should be proportional to (0| cos/3(/)|i?4) which 
is a finite number. The one-breather form factors of the 
type (0|e*'^'^|i?„(0)) = Fb„ can be computed from the 
residue of the soliton-antisoliton form factors at points 
9„ = j7r(l — nS,)), or, equivalently from the bootstrap 
(the definition of Fb„ can be found in [1^). For exam- 
ple, the breather B2 is a bound state of two breathers 
Bi. Using this bootstrap procedure we can find 



F\ 



.cxp(2/30)_ 



where 



gi}^/2coi(^) sin(7rO exp[/(-6')]e' 

cot(^)n:=icot2(M£) 



m 



dt sinh^ [t (1 - f )] sinh[i(e - 1)] 



t 



sinh[2i] cosh[t] sinh[t^] 



(45) 



(46) 



The saturation of the form factors implies that the ma- 
trix elements (0| cos{(3(j))\k) behave as 1/E{k) at large k 



ensuring the convergence of integrals in Eqs. ([50]) . (PTI) . 
and (|32|) and thus the validity of the universal scalings. In 
the next section using the conformal perturbation theory 
we will also show that there are no UV issues with these 
scaling relations, thus the lowest form factor expansion 
is indeed justified. We note that the saturation of the 
breather form factors can be also understood from the 
fact that cos(/?(^) corresponds to the trace of the stress- 
energy tensor, whereas breathers correspond to the field 
(f) itself. This saturation is thus expected for all breather 
contributions corresponding to exciting more than a sin- 
gle pair of breathers. In Fig. [1] wc show the dependence 
of the BiBi and B2B2 breather form factors as a func- 
tion of the rapidity. Both plots indicate clear saturation 
with at the value corresponding to the square of the 
single breather form factor. 



V. ANALYSIS OF THE HIGH ENERGY 
EXCITATIONS BASED ON THE FIDELITY 
SUSCEPTIBILITY AND THE CONFORMAL 
PERTURBATION THEORY 

The other approach for finding the scaling of Pcx and 
other quantities with the quench amplitude A/ is based 
on the analysis of the fidelity susceptibility x/j which is 
in turn related to the correlation functions. From Eq. ^ 
we know that for sudden quenches of small amplitude we 
can write: 



Pex(A/)«L|A/px/(A/), 



(47) 
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FIG. 1: BiB\ (top) and B2B2 (bottom) breather form factor 
for f = 1/41 as a function of the rapidity 9. The saturation 
lines are given by (Jsi)^ and (-^32)^ respectively. 



with Xf defined in Eq. A more accurate analysis 

based on the Cauchy-Schwartz inequality (see Sec. IVIip 
shows that Pox is bounded by the integral of Xf 
as Pcx(A/) < LA/ /p^ rfAx/(A), still the scaling of 
A/ x/(A)rfA and AyX/(A/) is the same for any power- 
law dependence of X/(A/), which we typically find near 
the critical point. Moreover in this section we will be 
primarily interested in the UV limit, corresponding to 
exciting high energy states, where Eq. (j47|) is justified 
by the conformal perturbation theory. Likewise one can 
obtain the expression for heat density: 



with 



Td H 



mv\n)? 



Po(A)' 



(48) 



(49) 



The susceptibility xe is equal to the generalized adia- 
batic susceptibility of order one: xe = Xi according to 
Eq. (0). 

For sudden quenches the wave function does not 
change during the quench. Hence the total energy af- 
ter the quench is given by the expectation value of the 
new Hamiltonian taken with the initial wavefunction: 



(50) 



where Hq is the free conformal Hamiltonian correspond- 
ing to the critical point and V = —2 J dx cos(/30). In the 
conformal limit the ground state is invariant under the 
uniform phase rotation (j){x) — >■ (l){x) + (f>o., where (f>o is an 
arbitrary phase. Using this symmetry the second term 



in Eq. ((50)) drops out and thus the total energy does 
not change during the quench. The heat is defined as 
the difference between the total energy and the adiabatic 
(ground state) energy of the instantaneous Hamiltonian 
{Egs{Xf)). Therefore we see that in our case the heat is 
equal by magnitude and opposite by sign to the change of 
the ground state energy: Q = — (i?gs(A/) — i?gs(0)). One 
can recognize that Eq. pS)) indeed gives (up to the sign) 
the second order perturbative correction to the ground 
state energy. As in the case of Pex, the perturbative ex- 
pression is justified if we are interested in understand- 
ing the UV limit, where the energy denominator (see 
Eq. (gg])) is large. 

As we argued in Sec.|ll]the susceptibilities xe and Xf 
can be expressed through the integral of the connected 
part of the imaginary time correlation function G{x,t): 



G{x,t) = (O|cos(/30(a;,T))cos(/30(O,O))|O) 
- (0|cos(/3(/.(0,0))|0)2. 



(51) 



Namely, for a translationally invariant system we have: 



XE 

Xf 



dr I dx G{x,t), 



I dr I dxTG{x,T). 
Jo Jo 



(52) 
(53) 



Let us point out that G{x, t) is a monotonically decreas- 
ing non-negative function of r = v^^ + r^ vanishing in 
the limit r — > 00. The non-negativity of J dx G{x,t) 
follows from the Lehmann's representation: 



/ dxG{x,T) ^ yY^mVln)]^ exp[-TiE„ - Eo 



)]• 

(54) 

The relativistic invariance of the SG model also implies 
the non-negativity of G{x, r) itself. 

The scaling of Pcx and Q with A/, at small A/, is 
thus related to the small A/ behavior of the integrals of 
G{x,t). The long distance behavior of the correlations 
functions at scales r 3> is given by the form factor 

expansion and can be analyzed by the methods of the pre- 
vious section, leading to essentially the same predictions. 
In particular, because at any finite A/ the spectrum of the 
SG model is gapped, G{x, t) exponentially decays on a 
correlation length scales r ss l/M^ ^ l/IA/p"*"^^/^. 
The scaling dimension of the operator cos{P(j){x,T)), en- 
tering the correlation function, is 0^ /Att — 2^/(1 -I- ^), 
implying that the scaling dimension of G is twice that. 
Thus we can estimate the contribution to the scaling 
of Xf coming from r, x ^ C, as Xf ^ ^3^(^^44/(1+^) ^ 
|A/|'-^+'S)/^~^ ^ |A/|''~^. This is precisely the scahng we 
obtained from the form factors expansion. 

The primary goal of this section is to analyze the UV 
(short distance) limit where the form factor expansion in 
the lowest order does not reproduce the correct scaling 
behavior, as we pointed out in Sec. IIVI For the contribu- 
tion to the susceptibilities and thus to Pcx and Q (as well 
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as to nex and Sd) coming from r <^ ( ^ one can use 
the cosine term as a small perturbation over the free (or 
conformal) limit. We note that by the uncertainty princi- 
ple the short distances correspond to exciting states with 
high energies (and momenta) En ^ Mg (fc ^ 1/C). In 
turn the short distance behavior of the correlation func- 
tions is governed by power-law asymptotics with the ex- 
ponents given by the scaling (conformal) dimensions of 
the operator cos{f3(j)), (which is 0^ /An = 2^/(^ -I- 1) as 
previously said). Thus in Euclidean notations; 



G{x,t) 



1 



(a;2 + r2)/3V4x■ 



(55) 



While this result is natural, it is somewhat nontrivial to 
obtain from the form factor perturbation theory since it 
requires an infinite re-summation of the form factors as 
it is nicely demonstrated in Refs. [1, UHl- So we can con- 
clude that in order to correctly reproduce the high energy 
contribution to Pqx and Q it is not sufficient to restrict 
the excited states in Eqs. (l28t and (|49)) to simply 
pairs of solitons and antisolitons with opposite momenta. 
It appears that only in the free fermionic point ^ = 1 
the lowest soliton-antisoliton contribution to G{x^ r) cor- 
rectly predicts the exponent in Eq. (|55|) . We believe that 
in general in our out-of-equilibrium setup it is possible to 
use a combination of some RG-type computations with 
the form factors based approach similar to the one in 
Ref. j66i]. 

From Eqs. ([5^ and we see that at the criti- 

cal point xe has IR divergence at /J^ < 47r correspond- 
ing to ^ < 1 and X/ has such divergence at /J^ < 67r 
corresponding to ^ < 3. By now we understand that 
this divergence simply implies that the scaling of Q{Xf) 
and Pex(A/) is non- analytic in A/ or the system size in 
these domains. This divergence also indicates that the 
main contribution to the susceptibilities, and thus to the 
heat and the probability of exciting the system, comes 
from small energies E ^ Mg, where one can reliably use 
the lowest form factors expansion. Conversely for ^ > 1 
the susceptibility xe has a UV divergence, which is cut- 
off by the non-universal short distance regularization of 
the SG model. Similarly for ^ > 3 the fidelity suscepti- 
bility X/ is dominated by short distances or high ener- 
gies. This implies that in these regimes of UV divergence, 
XeW and X/(A) approach non-universal constant values 
as A — >■ (note that these values are strictly positive be- 
cause G{x,t) is strictly positive). Therefore the leading 
asymptotical behavior for the heat becomes quadratic for 
^ > 1: Q ~ |A/|2x_b(0). Likewise the scahng for P^x be- 
comes extensive for ^ > 3: Pox ~ LXjxf{0)- This exactly 
agrees with the crossover to the quadratic scaling antic- 
ipated in the Sec. |lll that occurs for {d + z)v > 2 for 
Q and for dv > 2 for Pcx ■ This crossover is directly 
analogous to the crossover between non-analytic and an- 
alytic regimes of scaling for the linear quenches in higher 
dimensions [s^. Similar crossover between perturbative 
and non-perturbative scalings at ^ = 3 occurs for the 
density of excitations and the entropy. 



While the short distance behavior of the correlation 
function given by the conformal limit (j55p is intuitively 
anticipated, finding corrections to this scaling due to the 
finite value of the coupling A requires using conformal 
perturbation theory [67l] . In the latter the starting point 
is the conformal field theory (free boson compactified on 
a circle of radius R in our case). The Hilbert space then is 
the Verma module - representation space of the Virasoro 
algebra. The different basis states in this approach are 
labeled by two integers (n, m) corresponding to the 
eigenvalues p± = n/Rzt m/2R of the momenta of the 
left/right movers ip{z),(p{z), with ^{x,t) = ip{z) + (p{z). 
The perturbed theory is constructed by adding vertex 
operators V(„ „i)(z,z) = exp[i{p+ip{z) + p^Cp{z))] to the 
free theory. The scaling dimensions of the chiral parts of 
these vertex operators are A± = pj-/2 so that the scaling 
dimension h of V(^n,m} is h = A+ + A_. The sine-Gordon 
theory is considered then as a perturbation of this 
free theory by the vertex operators with n,m — {±1, 0}. 



HsG — HcFT + V, 

y = -2A/ / {V(i,a){z.z) + V(^ifl){z,z))dx 
Jo 



(56) 



$(a:, t) = VATr(t>{x, t), R = ^ Atk j . 



(57) 

The scaling dimension of V is then 0^ /An. Since the 
Hilbert space is built by the action of the vertex oper- 
ators on the vacuum state, |n, m) = V(„^m)(0,0)|wac), 
(higher-energy states are built by the application of pow- 
ers of the chiral components of the Kac-Moody currents 
to these states) the conformal perturbation theory is 
simply a perturbation theory in this space, and the ef- 
fect of the perturbation is an expansion in powers of the 
(MsL)2~'\ where h = /An is the scaling dimension of 
the perturbation (so this expansion is in integer powers 
of A/). In Ref. 62] such expansion to the third power of 
Xf was applied to the correlation function. In particular, 
it was found that 



(a;,r) -j/3</.(0,0) 



|0) 



1 



CX 



1 



(58) 

where r = y/x'^ + t'^ and C is a constant. In the repulsive 
regime Aw < 13'^ < Stt where this expansion was derived 
the second term clearly decays slower at large r than the 
first leading term indicating that the Taylor expansion in 
A / can break down. This is of course anticipated because 
cos{(3ip) is the relevant perturbation. Conversely at small 
r the first term corresponding to the conformal limit is 
dominant. 

If we use this expansion to compute xe according to 
Eq. (|52l) then as we discussed above for f3'^ < Air (or 
^ < 1) the first term gives IR divergence which should be 
fixed by using the proper long distance asymptotics of the 
correlation function. On the other hand for An < /3'^ < 6tt 
(1 < f < 3) the first term gives a finite constant contribu- 
tion to xe while the second term proportional to A^ gives 
a divergence which again has to be fixed by the proper 
long distance asymptotics. For 3 < ^ < 5 the second 
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term gives a convergent contribution indicating that xe 
has analytic expansion in A/ up to A ^ but the following 
term cx Aj is expected to give a non-analytic contribu- 
tion and so on. Thus from the conformal perturbation 
we can anticipate that the asymptotic expansion for the 
heat in powers of A has the following structure 

g= E Cr^f + B\\f\^+'+o{\\f\'+^), (59) 

0<n<[(l+C)/2] 

where C„ are some non-universal constants which depend 
on the short distance cutoff and i? is a universal constant 
determined by the low energy (long distance) contribu- 
tion to Xe (and thus to the heat). At the special points 
1-1-^ = 2n the constant B is expected to have logarithmic 
divergence. 

Similar analysis can be done for Xf using Eq. ([53]) and 
the expansion (f58| . Now we observe that the first con- 
formal contribution diverges at < Gtt (^ < 3). As 
before this divergence is cutoff by using the proper long 
distance asymptotics of G{x, r) and results in the non- 
analytic scahng of XfW- For Gtt < < Ttt (3 < ^ < 7) 
the first term in Eq. (j58p gives a constant contribution 
to X/: while the second term diverges and thus results in 
a non- analytic contribution. For 7 < ^ < 11 we see that 
both the first and the second terms in (|58|) give conver- 
gent contribution to x/i but we can anticipate that the 
third term gives a non-analytic correction and so on. We 
can therefore anticipate the scaling for the probability of 
exciting the system: 

(l+«)/4 

Pe. = i E C„Af -(-BA^L4/(i+«)+o(|A/pi4/(i+«)), 

n=l 

(60) 

where as before (?„ are non-universal constants and B is 
universal and also expected to have log-divergence at the 
special points 1 -f ^ = 4?i. While we can expect a similar 
structure of the expansion for rxcx and Sd^ it is harder to 
analyze these quantities in the UV regime because they 
do not have an explicit representation through the cor- 
relation functions. We simply note that the crossover to 
the leading quadratic scaling for these quantities when 
^ > 3 (with probable log corrections for the entropy) is 
expected on general grounds as a manifestation of the va- 
lidity of usual perturbation theory applicable to the high 
energy excitations, where the probability of the transition 
to be excited scales as |a„p '-^ Xj. Since the low energy 
excitations coming from creating soliton-antisoliton pairs 
give ncx, Sd ~ |A/P"'"^^/^ we expect that the high energy 
excitations will dominate precisely when (1 + £^)/2 > 2. 

One can actually analyze the energy of an arbitrary 
state obtained by acting with the operator on 
the vacuum j^*) = ^'|0), as a function of Xr using the 
nonlinear integration equation approach [681 l69j . The 
corresponding expansion in powers of A/ reads: 

^ = -^i+'-^ + Bl + ±C^m-^^-'^-^ (61) 

^ n— 1 



where / = MgL, B = — tan(7r^/2)/4 is the bulk energy in 
the thermodynamic limit, c is the central charge (equal 
to 1 in our case), C„(^') are non- universal constants. The 
bulk energy becomes infinite at the point where ^ is an 
odd integer. At this points there is a value of n for which 
n{2-h) = 1, i.e. 13^ = 87r(l - l/(2n)) and C„ oo. The 
infinite contributions from the bulk part and from the " C- 
part" cancel giving the logarithmic singularity propor- 
tional to l\og{l). These expectations are consistent with 
the findings of Refs. [6^, [g^. In particular, it was found 
that at the special points — 87r(l — l/2n), correspond- 
ing to ^ = 2n — 1, the expression for the ground state 
energy gives logarithmic singularities. Exactly at these 
points the short distance expansion method of Konik and 
LeClair ^6^ gives the divergency of the n-th order of the 
perturbative expansion. Noting that for sudden quenches 
the heat is related to the ground state energy and thus 
to Xe we see that the expansion (1611) is consistent with 
the asymptotical form ([5^ suggested above. 

Note that in the attractive regime < 47r, ^ < 1 the 
UV asymptotics of the correlation functions always give 
long distance divergent contribution to both xe and Xf 
indicating that the corresponding susceptibilities are de- 
termined by the IR region where the lowest form factors 
expansion is accurate. This confirms that the general 
scalings obtained in Ref. [i^ (see also Sec. and ver- 
ified in the previous section are justified in the whole 
breather region. 

Let us point out that a very similar analysis can be 
extended to slow quenches with arbitrary exponent r (see 
Eq. The validity of the perturbative quadratic and 
extensive scaling of Pcx is determined by the convergence 
of the adiabatic susceptibility of order 2r-|-2 (see Eqs. dU 
and dZ])). In particular, for linear quenches, the relevant 
susceptibility X4(A) is always IR divergent in the regime 
of our interest, < < 87r, implying that the scaling 
of Pex in the case of linear adiabatic quenches in one 
dimension remains non-analytic in the whole gapped 
regime of the SG-model and the form factors analysis in 
this range is always justified. Thus in order to analyze 
the crossover between analytic and nonanalytic regimes 
of scaling we need to consider the generalizations of the 
SG model to higher dimensions. This can be easily done 
in the two limits = 1 and ^ <SC 1 corresponding to free 
massive fermionic and bosonic theories respectively. We 
will discuss these generalizations and the crossovers in 
the next section. 



VI. EXACTLY SOLVABLE CASES: FREE 
MASSIVE BOSONS AND FERMIONS 

In this section we will focus on two limits of the SG 
model where the dynamics can be analyzed exactly with- 
out the need to rely on the adiabatic perturbation theory 
or other approximations. These limits have also the ad- 
vantage to have straightforward multi-dimensional gen- 
eralization. First we consider the limit of /3 ^ 1. In 
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this case the cosine term in the Hamiltonian (1221) can be 
expanded to the quadratic order in the phase (j) and the 
problem maps into the free scalar massive bosonic field 
theory. The quench process corresponds to a sudden in- 
crease of the mass starting from zero. The quench dy- 
namics in this regime was recently considered in Ref. [Tj 
from a different point of view, analyzing the correlation 
functions and the asymptotic steady state. The other 
solvable limit corresponds to /? = 2^ {£_ = 1). For 
this value of /3 the problem deals with hard-core bosons 
or non-interacting fermions. In the language of fermions 
the quench process corresponds to a sudden turn on of a 
commensurate periodic potential, which opens a gap at 
the fermi momentum. In the following two sections we 
will discuss these limits in detail, analyze the scalings of 
the density of quasiparticles, heat, and entropy and show 
how those results agree with the general scaling predic- 
tions of Sec. ini For these solvable cases in Sec. IVIDI we 
will be able also to extend the exact results to quenches 
at finite temperatures. 



A. Elementary derivation of the scaling relations 
using the adiabatic perturbation theory 

Before proceeding with the analysis of specific bosonic 
and fermionic models we show how the scaling relations 
and the crossover to the quadratic scaling can be under- 
stood if we are dealing with a free theory. This discussion 
goes along the lines with that of Ref. [29] for slow linear 
quenches. Using the general discussion of Sec. |IT]we find 
that within the adiabatic perturbation theory the am- 
plitude for creating an excited state in a sudden quench 
is: 



a„(A/) 



{n\V\Q) 



(62) 



If we are dealing with a free spatially uniform Hamilto- 
nian then the operator V must be quadratic in the cre- 
ation and annihilation operators and conserve the total 
momentum. Thus acting on the vacuum state V can ei- 
ther leave it intact or create a pair of quasiparticles with 
opposite momenta. Hence all the intermediate states can 
be characterized by the momentum label q. Thus the to- 
tal probability to excite the system is: 



Pcx = 1-11(1 

q>Q 



(63) 



where the product is taken over different pairs of two- 
particle states with opposite momenta and Pexiq) = 
|Q;g(A/)p is the probability to create the pair. In the 
limit when ^g^QPexil) ^ 1 the product can be approx- 
imated by the sum and we have: 



q>0 



(64) 



Note that since each term in the sum corresponds to a 
pair of quasiparticles the number of excited quasiparticles 



is given by the same expression above multiplied by a fac- 
tor of two, which can be absorbed into the sum by sum- 
ming over all positive and negative momentum states. 
Thus 



dX 



{q\V\0) 



Eq{\)-Eoi\) 



(65) 



where \q) is a short-hand notation for the state with two 
excited quasiparticles with momenta q and —q and Eq{X) 
is the corresponding excitation energy of the pair. Let us 
analyze the expression above using the scaling arguments 
(see also Refs. [1^ HO, Near the critical point we 

expect that Eq{\) - Eo{X) = \X\'"'F{q/X''), with F{x) a 
scaling function. In the limit a; 3> 1 we must have F(x) ~ 
to ensure the scale invariance at the critical point. 
Similarly we write the scaling of the matrix element: 



Eq{X) 



(66) 



where G{x) is another scaling function. The same re- 
quirement of scale invariance at the critical point gives 
that G{x) ~ 1/x^^'^ at a; ^ 1. Some examples where 
the scalings and asymptotics for the matrix element and 
the ener gy were explicitly confirmed can be found in 
Refs. [sol l55||. Using these scaling results it is straightfor- 
ward to find the scaling for the density of excited quasi- 
particles. In particular, focusing on the limit |A/| » 
L"^/"^, where the summation over momenta can be re- 
placed by the integration, and changing variables q — 
I A/ 1 "77 and X — Xf£^ in Eq. (|65|) we find that 



{2TrY 



s 



A similar expression is valid for the heat density: 



d'^r] 
(27r)° 



(67) 



(68) 



We thus see that we are reproducing the scaling predic- 
tion of Sec. ini Eqs. (|T3)) and (|T7)) . for sudden quenches 
(r = 0) as long as the limits of the integration over rj 
can be extended to infinity. This is only justified if the 
integral over ij converges at large 77. Using the asymptot- 
ical expression for the scaling functions G and F we see 
that this is the case for dv < 2 ii we are analyzing ricx 
and for {d + z)v < 2 if we are analyzing the heat. Alter- 
natively the integrals over the momenta are dominated 
by the high energy cutoff and the scalings of the corre- 
sponding quantities become quadratic with A/, in accord 
with the predictions of Sec. [Hi A similar analysis shows 
that the adiabatic perturbation theory reproduces the 
correct scaling also for the entropy density Sd ~ 
when dv < 2. The same approach confirms the correct 
scaling of all these quantities for small quench ampli- 

^cx,Pex ^ |A/|2L2/-d_ The 



tudes |A/| < l/i^/", e.g 
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easiest way to see this is from matching the quadratic 
scahng and the universal scaling when the correlations 
length ^ ~ 1/1 A/ 1" becomes of the order of the system 
size. In the next two sections we will explicitly show how 
these crossovers emerge from the two free limits of the 
SG model. 



Free massive bosons 



Consider the SG Hamiltonian (|22|) when /3 — ^ 0. In 
this limit we can Taylor expand the cosine term up to 
the quadratic order in the field 4> and obtain a quadratic 
Hamiltonian, that in the Fourier space has the form: 



H = 



Kq{\)\(t)q 



(69) 



with Kq{\) = + 2X{t)(3^ and Ht,0t = H.,,^-,. We 
point out that here as in Eq. (|22|) we have set hvs = 1, 
with Vs being the sound velocity. For each momentum 
q the fields HJ and (j)q are conjugates, [Hj,0q'] — i5q,q' . 
Therefore we are dealing with a sum of independent har- 
monic oscillators with a time-dependent compressibility 
K (or the inverse mass). The wave function for this time- 
dependent Hamiltonian can be found exactly, for instance 
in Refs. [501,^5^ this was done for the case of an adiabatic- 
linear quench (A(i) = vt). Here we focus on the case of 
a sudden quench and show that the exact results agree 
with the general scaling predictions that have been pre- 
sented in the SecHIl Furthermore the Hamiltonian ([M)) 
can be easily extended to higher dimensions and thus we 
are able to see how the crossover to the quadratic scalings 
and the usual perturbation theory emerges as the corre- 
sponding exponents dv or (d+ z)v exceed two. We note 
that in this limit of the SG model v = 1/2 and z = 1. 
For the free theory these exponents remain the same in 
all dimensions. 

From an experimental point of view this limit can be 
realized in the case of merging (coupling through tun- 
neling) two one-dimensional Bose gases; indeed as it has 
been pointed out in [l^ the SG hamiltonian (|22p de- 
scribes themerging process through the identification of 
/3 = yJOjijK, with K the Luttinger parameter for each 
Bose gas, and X{t) proportional to the tunneling strength 
between the two tubes. In the limit of weakly interacting 
gases K ^ oo, hence /3 — and the process is described 
by the Hamiltonian (IMl) . 

Let us consider a sudden quench in the Hamiltonian 
([69| where A is changed in time as: X{t) = XfQ{t), where 
Q{t) is the step function. It is convenient to switch to the 
representation in terms of the creation and annihilation 
operators a' and a, such that [a, a^] = 1, in order to work 
with the number-state vectors: |n) — {a^ / Vn\\0) . In 



particular, for a fixed A those operators are: 

""^^^^ = 77T^ + *\/^ '^0 , (70) 
(4k(A))4 V ^ / 

at(A)=(^(nt-.V^0t), (n) 

for each mode q, such that the Hamiltonian in this basis 
becomes: H = Eg V i^Wi^lWaqiX) + 1/2). For uni- 
form quenches all g-modes are independent from each 
other and the dynamics of the system is factorizable. 
We are interested in evaluating the effects of the sudden 
quench on the population of the excited levels, therefore 
we want to express the initial ground state |0g)o before 
the quench (A = 0) in terms of the eigenstates {\nq)xj,} 
of the final Hamiltonian (A = Xf): 



10, 



'^C2niq)\2nq)xf, 



(72) 



where we used the fact that only even modes are excited 
in the quench process because of the parity conservation. 
The coefficients C2n{q) can be evaluated expressing the 
operator ag(0) in terms of aJ(A/) and ag(A/): 



C2n{q) = (-1)"1 



'(2n-l)!! / v/MATT- VMO) 



(2n) 



VKg(A/) + y^70) 



2{nq{Xf)KqiO)y 



(73) 



One can check that these coefficients are properly nor- 
malized: En |c2n('z)P = 1- Physically each state \2nq) 
represents 2nq quasiparticles with opposite and equal mo- 
menta q and —q. Therefore the average number of exci- 
tations in each mode is nex{q) = En 2'^|c2n('?)P- We can 
also formally define the probability that the mode with 
momentum q is excited as: 



Pcx{q) = ^ |c2n(g)|' 

n=iQ 



(74) 



The total probability of exciting the system Pcx is found 
as the complementary of the joint probability that no 
mode is excited after the quench: 

Pcx = 1 - n M = 1 - 11(1 - pex(9)). (75) 

g>0 <j>0 

In the limit EgPcx('i') <C 1 we see that Pcx{q) ~ 
EgPcx(g)- 

From Eq. ([75)1 we find the condition Pcx('z) <C 1 is 
satisfied if we impose: 



7Ma7)~ VMo) ^ V¥ + 2XJ^~q 
v/M^+v/^o) + + q 



< 1 (76) 



that implies 2Xf0^ <C q^ ■ In this limit we have Pex{q) ~ 
|A/p/3^/(8(7''). For a system with size L the lowest 
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nonzero momentum g is ^ and therefore the probabih- 
ties that Pcxil) ^ 1 for all q is equivalent to: 



A/ < A, 



1 /27r 



(77) 



It is obvious that this condition also ensures that 
J2qPex{q) ^ 1- In this limit the total probability for 
the system to be excited is: 



8(27r) 



k(4), (78) 



where ^(fc) is the Riemann's Zeta function. It is easy to 
check that the coefficient multiplying iA^ is precisely the 
fidelity susceptibility Xf evaluated in the massless limit 
A = 0. The super linear scaling of the excitation proba- 
bility Pex IL ~ is consistent with the one anticipated 
in Ref. [H (see also SecHU: x/(0) ~ L-'^+^Z" = L^. 

For quenches of larger amplitude, when A/ > As, as we 
already discussed in Sec|ll]the probability of exciting the 
system Pex is almost unity. This indicates that Pex can 
not differentiate between different excited states. This 
comes from the fact that even if we excite a single pair 
of quasiparticles, the state becomes immediately orthog- 
onal to the ground state and Pex = 1. At the same time 
physically this state is almost indistinguishable from the 
ground state since a single pair of excited quasiparticles 
can not affect any thermodynamic properties of the sys- 
tem. Let us analyze instead that following object: 



EaPox(9) 



1 



0.036 



(79) 

which represents the sum of probabilities of excitations of 
different momentum modes. In general such object can 
be defined first evaluating the two-particle density ma- 
trix in the momentum space and then finding probabil- 
ity that the corresponding mode is excited. To evaluate 
Eq. ([7^1 we used the fact that the sum over momenta 
can be converted into the integral and for small ampli- 
tude quenches the limits of integration can be extended 
to [0, cx)). The scaling dependence in (|79p agrees with 
^ |A/Px/(A/) = IXfl'^" and hence it illustrates that the 
fidelity susceptibility actually describes the behavior of 
"^qPcxig) rather than Pex. Only for infinitesimally small 
quenches |A/| <C |As|, when these two objects coincide, 
the fidelity susceptibility also describes the scaling of Pex- 
This point can be made even more explicit if we consider 
a quench starting and ending in the massive limit Ai > 
and Xf = Xi + SX > 0. Then one can check that as 
long as \SX\ < A^ we have X]„Pcx(g) « L{SX)'^Xf{X)- 
At the same time a similar scaling for Pox is only valid 
for \6X\ < 1/VL i.e. in the vanishingly small interval of 
quench amplitudes in the thermodynamic limit. 

The sum of probabilities in ((7^ itself is not a phys- 
ical observable. A measurable quantity with the same 
scaling is the density of created quasiparticles nex = 



l/-^I]q'^cx(g), where n^xiq) = {a\aq). In the limit 
(g) ^ 1 the sum in ([71)) is dominated by the first 



P' 

term, n = 1, corresponding to the creation of one pair 
of quasiparticles with opposite momenta. Thus in this 
limit ricxiq) ~ 2pex(9)- From this we conclude that for 
the small amplitude quenches A ^ Ag we have riex ~ 2pcx, 
with = Pcx/L given by Eq. (|75)) . As we will see the 
same relation is true if we are dealing with fermionic sys- 
tems. Physically this relation comes from the fact that 
we are dealing with few-body operators. Thus in the low- 
est order of perturbation theory the coupling through the 
quench dynamics leads to generating a small number of 
quasiparticles in each momentum mode. If the system 
is noninteracting (or more generally integrable) then the 
number of these quasiparticles is conserved in time af- 
ter the quench. In nonintegrable systems nex changes in 
time due to scattering processes. Therefore the scaling of 
ricx with A/ is expected to be valid only either at times 
shorter than the relaxation time, or if the generated de- 
fects are topologically protected like in the Kibble-Zurek 
mechanism [s^, IS] • We note that the scaling of the heat 
is insensitive to the relaxation processes because the en- 
ergy in a closed system is conserved after the quench 
process. The same is true about the (diagonal) entropy. 

For finite quench amplitudes A > A^ the expressions 
for riQxiq) and ^/ L^^pcx{q) are different. In particular, 
it is easy to check that 



(V>%(a7)- Vn(o))' 

4v/«,(A/)«,(0) 



(80) 



For q > 2A/;32 ^hat n^xiq) ~ \X\^I3'^ /{iq^) « 

2pcxiq)- In the opposite limit Uexiq) ~ \J'2XpP- jq while 
Pcxiq) ~ 1- The difference between ricxiq) and Pcx(<z) 
happens because the states with more than two quasi- 
particles per each mode are excited for small momenta. 
Such states contribute differently into the probability of 
excitation and the quasiparticle density. For the total 
density of quasiparticles we thus find: 



^y^dgnex(g)«— ^log 



A^L2 
27r2 



(81) 

where A is the upper momentum cutoff. Note that this 
integral gives additional logarithmic dependence on both 
the quench amplitude and the system size. This depen- 
dence appears because of the infrared divergence of the 
density of excited quasiparticles at small q. It is very sim- 
ilar in nature to the onset of the non-adiabatic regime for 
slow quenches [s^ and is due to the bunching effect. As 
we will argue later this divergence disappears in dimen- 
sions higher than one where the scalings of both riex and 
TliqPcxiq) I L''' with the quench amplitude become identi- 
cal. In one dimension the two scalings agree apart from 
a weak logarithmic correction. 

Similarly we can analyze the heat density, or equiv- 
alently the non-adiabatic energy change in the system 
due to the quench. For A < A^ we can find, similarly to 
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Eq. (EID: 



1 i2^4^2 

L ^ 

g>0 g>0 



L ^ 4q3 4(27r)3 



<(3). 

(82) 



In the opposite limit A > As we have: 
2t 



1- r + 2A//32ne.(g) 



27r 



27r2 



(83) 



Up to the logarithmic correction in the above equa- 
tion, these results agree with the prediction in Sec. |lll 
Q - >^^Xe{X) with xe{X) ~ L2/i.-(d+^) foj. A < As and 
Xb(A) ^ \\\id+^>--2 for A > As. We point out that if 
we consider an opposite process where we start with a 
small initial coupling Ai and end at the critical point, 
then the expression for will remain the same, while 
the expression for the heat for A > As will change because 
the energy is now evaluated at the critical point. This 
will remove the extra logarithmic divergence and we will 
simply get Q oc Ai/?^. 

Finally let us consider the density of the (d-)entropy 
generated in the quench. Like Pex it gives the measure 
of the excitation of the system, however, it is explicitly 
extensive and has well defined scaling properties both 
for A < As and A > As. Unlike Uox, the entropy (like the 
heat) can be defined for any system, integrable or not. 
Because all momentum modes are independent we find 
that Sd = l/iZ]g>o where 



= -X|c2„(g)|'log(|c2„((?)|^ 



(84) 



For A ^ As the entropy of each mode is dominated 
by the lowest excitation: Sq « — 1 02(9) P log |c2((;)P| — 
-Pex(g)logPex(g), therefore: 



8(2^)4 



C(4) log 



8(2^)4 



(85) 



We note that for A < As the argument of the logarithm 
is smaller then one so the entropy is positive as it should. 
In the opposite limit A > As we find: 



Sd « 0.14 



(86) 



So we sec that the entropy has the same scaling as riox 
and l/LY,qPc^{q)- 

A similar analysis can be done also in the case of a 
linear quench, r = 1 and v = 5. As it has been shown in 
Ref. [50| , the scalings for the density of excited quasipar- 
ticles: nex ^ 5^^^ and for the heat, Q ~ 6"^^^ , in the case 
of quenches with rate 5 ^ l/Li^ agree with the predic- 
tions in Eqs. and (IT7|) . Similarly it is easy to check 
the correct scalings in the other regime, 5 <^ \jL^. For 



instance in this regime using the expression for the num- 
ber of excitation per mode q (see Eq. (103) in Ref. 50] ), 
we obtain: 



~mf ~ 16(27r)6 



C(6), (87) 



1/2, 



confirming the general scaling (1121) . which for v 
d— 1, z — 1, and r —\ gives Hex ~ L^S^ . 

Generalization to higher dimensions. Let us note that 
the Hamiltonian (I69|) can be analyzed in any spatial di- 
mension. The quench process simply describes the re- 
sponse of the free bosonic theory to a sudden turn on 
of the mass term. Since all the momentum modes are 
independent from each other the expressions for C2n{<l), 
Pcxil), 'T'cx('z) and Sq remain the same. The only dif- 
ference with one dimension is that we have to sum over 
all modes in the d-dimensional space. This introduces an 
additional density of states factor p{q) oc q'^~^ into all the 
expressions. Let us focus only on quenches with A > As. 
It is easy to check that for 1 < d < 4 one can rescale 
the momentum q — ^2\jfP-r\ and set the limits of inte- 
gration over r\ to [0,oo) when analyzing 1/i X]gPcx('?), 
riox, and Sd- This follows from the fact that at large 
q we have Pex{q) ^ l/?^ so the integrals over momenta 
converge and the upper cutoff can be sent to infinity. 
Likewise in all dimensions larger than one there are no 
infrared divergencies in all the integrals so the lower limit 
of integration over q: qmin ~ 2Tr/L can be sent to zero. 
Then for all these quantities we get the same desired 
scaling, e. g. rtcx ~ Cd{\fP'^Y^'^ , where Cd is a number 
which depends only on the dimensionality. Similarly for 
the heat for d < 3 we get Q w C'd(A//32)('^+i)/2. Above 
four dimensions (three dimensions for the heat) the inte- 
grals over momentum become ultraviolet divergent. So 
for c? > 3 the leading contribution to heat will be deter- 
mined by the large momentum asymptotic of the transi- 
tion probability: 



d"^ 



q q 



(88) 



where A ^ tt is the high-momentum cutoff and C is 
a non-universal constant. A similar expression is valid 
for ricx ~ ^/I^'^qPexiq) and for the entropy above 
four dimensions with the cutoff entering in the power 
Pqj. -(-jj^g entropy there is an additional loga- 
rithmic factor log(Aj/3''/A'*). In three (four) dimen- 
sions heat (density of quasiparticles) acquire an addi- 
tional logarithmic dependence on the quench amplitude: 
Q ^ Aj/?** log(A/(A//32)). This divergence indicates the 
crossover from quadratic to universal leading asymp- 
totics. In Sec. |V]we already discussed how similar loga- 
rithmic corrections show up in heat in one dimensional 
SG model when K = 1 for Q {K = 3/2 for riex)- 

Using the free bosonic theory it is also easy to verify 
that for 3 < d < 5 the first subleading asymptotic in the 
expansion of Q in powers of A/ is non-analytic and cutoff 
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independent. This follows from the fact that: 



Q = 



(89) 



Since nox(g) - X}P'^/4:q'^ \)I3^ I2q^ the first integral 

in the above equation is ultraviolet convergent below five 
dimensions and we can use the rescaling q — yJIXffS^ri 
and send limits of integration over ry to [0,oo). We can 
continue and get the asymptotical expression for the heat 
in high dimensions similar to Eq. ()59p 



Q 



E 

2<Ti<[(d+l)/2] 



(90) 

We note that this analytic expansion goes in both even 
and odd powers of A/ (unlike Eq. ([5^ ) where only even 
powers of A appear. This has to do with the fact that 
for the massive bosonic theory (unlike the SG model) 
A ^ — A is not a symmetry of the Hamiltonian. More- 
over, the theory is well defined only for A > 0. As in 
Eq. (|59p we find logarithmic corrections when nonana- 
lytic and analytic powers coincide, i.e. when (d + l)/2 
becomes an integer. A very similar expression holds for 
ncx with the only difference that in Eq. ((90)) one needs 
to replace d + 1 with d. We point out that for small 
amplitude quenches A < A^ the expansion (|90p will re- 
main valid with the first non-universal analytic terms 
being unaffected, while in the last non-analytic terms 
one needs to substitute |A/|(''+i)/^ to XjL^''^. There- 
fore above three dimensions for quenches of small am- 
plitude the non-analytic correction to the heat becomes 
sub-extensive. The same is true for the density of quasi- 
particles where the non-analytic term scales as XjL'^^'^. 



C. Free massive fermions 

Another important application of the SG Hamilto- 
nian (|22|) to physical systems feasible in cold atoms ex- 
periments, is given by the description of the loading pro- 
cess of a one-dimensional Bose gas into a commensurate 
optical lattice. In this case (3 is related to the Lut- 
tinger parameter through: (3 = I^/txK^ and the ampli- 
tude of the optical lattice V(f) is directly proportional 
to A [1, [13] • This process can be studied exactly in the 
Tonks-Girardeau regime, i.e. when K = \. In this case 
the repulsive interaction between the bosons is infinitely 
strong and the particles behave as impenetrable spheres 
(hard-core bosons) . It is well known that in this limit the 
system can be mapped into an equivalent system of free 
spinless fermions |70[. Therefore the dynamical problem 
of loading hard-core bosons into a commensurate optical 
lattice can instead be approached with the much sim- 
pler analysis of free massive fermions in a periodic po- 



tential. We note that this fermionic limit also describes 
the transverse field Ising model via the Jordan- Wigner 
transformation Q. Then the quench dynamics is per- 
formed by suddenly change the transverse field starting 
at the critical point. 

To understand the dynamics in the TG-gas limit, we 
need to solve the Schrodinger equation describing free 
fermions in a periodic commensurate potential with time 
dependent amplitude V{x,t) = T^(t) cos(2g/a;), where 
qj = IT /a is the Fermi momentum and a is the lattice 
spacing, which is also equal to the inverse particle den- 
sity. In the TG-limit the sound velocity Vs — Vf/K 
is equal to the Fermi velocity Vf. It is convenient to 
work in units such that a — 1, Vf = 1 and h = 1. The 
free fermion Hamiltonian becomes equivalent to the SG 
Hamiltonian if we additionally restrict the analysis 
to the two lowest bands of the Brillouin zone and lin- 
earize the spectrum close to the Fermi momentum. Then 
for each momentum the loa ding process is described by a 
Landau-Zener Hamiltonian (Tlj (see also Ref. [s^l for the 
discussion of the equivalent transverse field Ising model) : 



V{t)/2 k 
k -V{t)/2 



(91) 



where we defined the momentum k = q — qf, so that 
is measured from the Fermi momentum. By matching 
this spectrum with the one of the SG model we find that 
A = V/{Att). The extra factor of An appears because of 
the details of the bosonization procedure [t^, [t^ ■ 

The Hamiltonian ([?!]) can be diagonalized on the ba- 
sis of its eigenvectors: |-|-fc) and \—k), with respectively 
positive and negative energy. To study the probability of 
exciting the system due to a sudden quench of the cou- 
pling V: V{t) = Vf9{t), we need to evaluate the overlap 
of the initial ground state |— fe)o with the final excited 



state 



-k}Vr 



Therefore Pox(fc) 



ticular we find: 



-fe|-fe)oP- In par- 



Pox (fc) = 2 



+ fc2 _ Vf\k\ + {Vf - |fc|)^V7 + fc2 



V^ + k^ + VfJV^ + k^ 



(92) 

This transition probability has the following asymptotics: 
Pcx(fc) « V}V(4fc2) for fc » y/ and Pex(fc) « 1/2 for 
k ^ Vf. As in the previous case of free bosons, the ele- 
mentary excitation correspond to the creation of particle- 
hole pairs moving with opposite momenta (recall that 
annihilating a particle in the lower band can be viewed 
as creating a hole there with an opposite momentum). 
However, unlike with bosons, we can not create more 
than one pair for each momentum state due to the Pauli 
exclusion principle. This means that l/iX]fcPcx(fc) and 
ricx are always identical up to a factor of two. As in 
the bosonic case the crossover scale Vs can be defined 
from the requirement that X]fc^'cx(fc) ^ 1, which gives 
Vs ^ A-k/L (or equivalently A^ ^ l/£). For V < Vs 
the total probability to excite the system is small so that 
Pox — X]fc?'cx(fc)- This characteristic scale agrees with 
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the expected one As ~ L~^/'^ from the analysis of the 
system close to the quantum critical point, since in this 
case v = 1. Using the explicit asymptotics for the exci- 
tation probability Pcx{k) at small amplitudes V/ we find 
that in this case: 



3 f 



L 2 48 
For large quenches Vf ^ Vg we find that 



° dk 



Pcx(fc) 



4A/. 



(93) 



(94) 



We see that as in the case of bosons Uox has the correct 
scaling nex ~ |A/|'''' for A > and nox A^L^/''"'' for 
A < As in agreement with Eqs. and ((T5)l respectively. 
Likewise one can evaluate the heat density Q = 

2/Lj2kPeAk)Jk'^ + V7/4. For A < As wc find: 



-^y'fc2 + y//4pex(fc) 

k 

^log{L/2) = 8TTX}\og{L/2), 



(95) 



and in the opposite limit A > As 
Q^^\og{2n/Vf)^- 



-27rAjlog(2A/). (96) 



This scaling of heat agrees with the general expectations 
of Eqs. ([T5| and P7|) and the results obtained in the 
previous sections using the adiabatic perturbation theory. 
Since we have {d+z)^ = 2, there is an additional logarith- 
mic correction indicating the crossover from quadratic to 
nonanalytic scaling in the leading asymptotic of the scal- 
ing of Q with A/. Note that the difference in scaling 
between A < As and A > As only appears in the logarith- 
mic term. 

Finally we can evaluate the entropy density. Let us 
explicitly quote only the result for A > As: 



Sd = -^X![[^'='^(^)^°s(Pcx(fc)) 



+ (1 - Pex(fc) log(l - Pex(fc))] = Y - 27rXf. (97) 

The scaling of the entropy is thus again in accord with 
the exact solution. It is easy to check that for A < As 
the entropy density scales as S* ^ A^Llog(A/i), i.e. it is 
super-extensive. 

As in the bosonic case, these predictions can be also 
checked in the case of slow quenches. In Ref. [sOl it has 
been shown that for the linear quench X{t) = St one finds 
JT-cx ~ VS for S >• l/L^, in agreement with Eq. ([T3l) . 
while for S ^ one can find that (see Eq. (108) and 
Appendix A in [soj for the expression for Uq): 



k 



IE 



452 
1^ 



(2^)4^^ 



(98) 



which confirms the scaling ricx ^ L^5^ given by Eq. (|12l) . 

Generalization to higher dimensions. Like in the case 
of bosons one can extend the analysis to higher dimen- 
sional free fermionic models. For example, the two- 
dimensional generalization of the free fermionic theory 
describes commensurate fermions in the honeycomb lat- 
tice, which are realized in graphene [t^]. Such lattices can 
be also in principle realized in cold atoms using multiple 
laser beams [tI] . The quench process corresponds to the 
distortion of the lattice, which opens a gap in the spec- 
trum. Quench dynamics in similar setup was recently 
analyzed in Ref. [t^ . 

Unlike in the bosonic limit, where i' — 1/2, here we 
have v ^ 1. This implies that we expect that the 
quadratic scaling for the quasiparticle density and en- 
tropy emerges above two dimensions and for the heat 
above one dimension. This is indeed the case since 
the excitation probability at large momenta scales as 
Pcx{k) ~ l/fc2 (versus l/Zc^ in the case of bosons). To un- 
derstand the general structure of the heat density above 
one dimension we can more closely examine the expan- 
sion of the transition probability Pcx{k) at large k. In 
particular, from Eq. (j92p we find: 



Pcx(fc) 



1/2 



Wf 



(99) 



4fc2 I6fc4 32fc6 

As we see the expansion is in even powers of Vf which 
reflects the symmetry of the problem with respect to 
Vf — —Vf. For 1 < d < 3 we can rewrite the heat 
as 



Q « 2M 



2M 



° d^k V] 
_^ {2tiY 4fc ' 



^//4 ( Pex(fc) - ^ 



(100) 



where M is the number of independent light cones in the 
system {M — 2 for graphene). The second contribution 
gives a non-universal quadratic term, while in the first in- 
tegral we can rescale the variables fc — V/r/ and send the 
limits of integration over r\ to [0,oo). This contribution 
thus gives the non-analytic term jV/l'''''"^^'^. A similar 
analysis can be extended to higher dimensions and as a 
result we get the following expansion: 



Q 



E 

l<ri<[(d-|-l)/2] 



C„|A/|2"-|-i3|A/|'^+i-fo(|A|'^+i). (101) 



The coefficients C„ are non-universal (cutoff dependent) 
while the coefficient B is universal and can be obtained 
from the low energy description of the system. Similar 
considerations can be applied to the density of quasipar- 
ticles resulting in the expression similar to Eq. (llOip with 
rf + 1 — > d. For small amplitude quenches A < As only 
the universal term in the expansion above gets modified 
with |A/|'^'^^ — T' X'jL^~'^. As we already mentioned in the 
points where the power of the non-analytic term crosses 
the one of the analytic contribution {d + 1 = 2n) there 
are additional logarithmic corrections to the heat. 
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D. Finite temperature quenches 

Let us extend the analysis of these thermodynamic 
quantities to the case of a finite initial temperatures. 
We suppose to perform the same quench process but 
now preparing the system initially in a thermal state at 
temperature T. We will assume that the system is iso- 
lated during the dynamical process. For instantaneous 
quenches this assumption is generally always satisfied. 
For slow quenches this assumption is only valid if we are 
interested in the time scales shorter than the relaxation 
time with the thermal bath. This is the usual require- 
ment in order to be in the adiabatic limit in conventional 
thermodynamics [77j and it is usually well satisfied in 
cold atoms. We note that other studies have addressed 
the case of systems where the contact with a thermal 
bath during the dynamical process is essential [zl]. The 
advantage of focusing on the setup where the tempera- 
ture enters only through the initial conditions and does 
not affect the equations of motion is that one can expect 
that the scaling of various quantities will remain univer- 
sal. We note that finite temperature dynamics in this 
setup were studied earlier in Refs. [13, H^l for slow linear 
quenches and in Ref. [l^ for sudden quenches. 



1. Free massive bosons. 

For the case in consideration of a system of indepen- 
dent harmonic oscillators (see Eq. a finite initial 
temperature leads to the occupation of the eigen-modes 
according to the Bose-Einstein distribution: 



1 



1 



exp[/3V'^;(0)] - 1 exp[/?g] - 1 ' 



(102) 



This initial occupation enhances the transition probabil- 
ities. In Ref. [55.] it was shown that for an arbitrary time 
dependence of X{t), the mode occupation with momen- 
tum q at finite temperature is related to the one at zero 
temperature via: 



'(g) = ncx('7)coth 



nox((7) coth 



v/^O) 



2T 



nriq) 



2T 



nriq). 



(103) 



where n*°' {q) is the total number of excited quasiparti- 
cles during the dynamical process, that is indeed equal 
to the initial thermal population of this mode (71^(9)), 
plus the number of particles generated in the identical 
process at zero temperature multiplied by a hyperbolic 
cotangent factor. Therefore we can define the excita- 
tions n^(g) created only by the dynamical process at a 
finite temperature as the difference nl'^{q) — nT{q). It 
is immediately clear that the hyperbolic factor gives an 
enhancement of the transitions since at T 3> g it scales as 
2T/q ^ 1. The result (|103p can be seen either from the 



thermal average of the energy of each mode q as it was 
recently pointed out in Ref. [l^, or writing the initial 
thermal state in the Wigner form, as shown in Ref. [sH ). 

Let us now analyze how this additional factor affects 
the scaling of the quantities of interest. To simplify the 
discussion we will focus only on the heat density. We will 
also assume that the temperature is large compared to 
both the energy associated with the quench amplitude 
T y/\f(3'^ and the finite size quantization T 1/L. 
Otherwise we will return to the zero temperature asymp- 
totics. 

For A < As by simple modification of Eq. we find 



Q 



2T ^ ^^X} 



q>0 



4q4 



2(27r)' 



^C(4). 



(104) 



In the opposite limit A > As instead of Eq. we find: 



2r^ 



2Xf/3' 



q>Q 



-TLXff3'. 



(105) 



For small quenches the result (|104l) is fully consistent 
with the scaling prediction of the Eq. ([TB)) : Q ^ 
A^L^/"^"^. While for larger quench amplitudes we get 

a deviation. Instead of Q ^ \Xf\'^'^ = \/\Xf \ we get the 
scaling Q ~ |A/|L. The origin of this discrepancy is 
the infrared divergence of the sum in Eq. (|105p at small 
momenta. This divergence is similar to the one appear- 
ing in linear quenches, which leads to the non-adiabatic 
regime [55)]. As it is easy to check this divergence dis- 
appears in higher dimensions and the predicted general 
scaling is restored. Indeed for 2 < d < 4 and for A > As 
one can substitute the summation over momenta in the 
general expression for q by an integration, change the 
variables q = ^j2XfP'^ri and extend the limits of inte- 
gration over 77 to [0,00). This immediately implies that 
Q ~ T|A//3^|'*/^. Above four dimensions the integral 
over momenta becomes ultra-violet divergent. In this 
case the transitions are dominated by high momenta for 
which the initial thermal occupation is not important (if 
the temperature is much smaller than the high-energy 
cutoff: T <C TT in our units) and we are back to the zero 
temperature non- universal result ()88|) . It is easy to check 
that in this case (as at zero temperatures) the universal 
nonanalytic correction to Q becomes subleading thus in- 
stead of Eq. ^ we find: 



Q 



2<ri<[d/2] 



C„|A/|" + BT\Xff^ + o{\Xf^). (106) 



To summarize the finite temperature results for the 
bosonic systems we see that the substitution d ^ d — z 
always predicts the correct finite temperature scaling of 
the heat for A < As. In the opposite limit A > As this 
substitution d — > d — z also works in two dimensions and 
above but fails in one dimension. The same applies to 
the density of generated quasiparticles. 
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2. Free massive fermions. 

It is also straightforward to obtain the finite temper- 
ature asymptotics in the fermionic case. Since in this 
case we are deaUng with a sum of independent two-level 
systems, at finite temperature each upper and lower 
|— fc) level is occupied according to the Fermi distribution: 



cxp 



T 



1 



where we used the convention that the Fermi velocity is 
one so that the energy of the two states e±{k) ~ ztk. It 
is well known that for the two-level system the transition 
probability satisfies the detailed balance, i.e. the transi- 
tion probability from the top to the bottom level is the 
same as for the opposite process (for systems with more 
than two levels the detailed balance is generally not true, 
see e.g. Ref. Therefore the number of additionally 

excited particles in the dynamical process gets corrected 
as: 



nJx(fc) = ncx(fc)(4 - fk) ^ ncx(fc)tanh 



2T 



(107) 



This expression is similar to the one for the bosonic case 
P03p . with now the hyperbolic tangent factor replacing 
the cotangent one. Contrary to the bosonic case there- 
fore, the effect of the finite temperatures is to suppress 
the transitions since I tanh(a;) I < 1. At high temperatures 
T ^ k the suppression factor is approximately |fc|/(2T). 

As in the case of bosons let us focus on the scaling of 
heat. Note that the main contribution to Eqs. ((95)) and 
(|96p even at zero temperature comes from high energies 
k ^ A. Therefore in the leading order in A/ the temper- 
ature will not significantly modify the expression for the 
heat except changing the cutoff in the logarithm: 



Q«87rA2log(VT). 



(108) 



This scaling is valid both for A < A^ and A > A^ . This is 
indeed an anticipated result because the universal power 
(d -I- 2z)j/ = 3 is bigger than two. The next subleading 
correction in A/ remains universal though and scales as 
Xj/L for A < As and |A/p/T for A > Ag. It is straightfor- 
ward to check that this general structure of the expansion 
extends to higher dimensions so that for A > A.; instead 
of Eq. (HnH) we find: 



Cn\Xf\''"+B\Xff+^+oi\X\''+''). (109) 



l<n<[(d-|-2)/2] 



The non-universal coefficients C„ are temperature in- 
dependent (for r ^ A) while the universal coefficient 
B is different from the zero-temperature one. We thus 
see that the finite temperature effects result in changing 
d ^ d + z in the universal part of the expression for heat. 
It is easy to check that the same is true for the density 
of excited quasiparticles. 



VII. QUANTUM GEOMETRIC TENSORS AND 
FIDELITY SUSCEPTIBILITY. MAZUR 
INEQUALITIES. 

A. Quantum geometric tensors: generalizations 

In the discussion above we only considered situations 
where A was a single component coupling. In principle 
one can extend the analysis to more general setups where 
A is a M-component vector in the parametric space. 
Within the adiabatic perturbation theory it is easy to 
generalize the expression for the transition probability to 
the many-body state |n) due to a sudden quench in the 
multicomponent case: 



l««(A)|^ 



dX'- 



{n\dx'H\0) 



E^{\') - Eo{\') 



(110) 



where the contour integral is taken over an arbitrary path 
connecting the initial and final couplings. For sudden 
quenches the precise path is clearly unimportant. Next 
we will use the Cauchy-Schwarz inequality, which for ar- 
bitrary square integrable functions f{x), g{x) states that 



f(x)g{x)dx 



< 



\fix)\'dx- / Igix^dx. (Ill) 



Choosing g{x) = 1 and f{x) to be the integrand in 
Eq. (fTTU)) we find: 

Pcx = J2 \^n{X)\^ < A„ / dX'^Qa.A>^'), 

n^O a=l •'0 

(112) 

where Qa.aiX) are the diagonal components of the quan- 
tum geometric tensor Qa.fiiX) defined in Ref. [H^. This 
tensor and its rescaled counterpart {qap = L~'^Qaj3) has 
definite scaling properties close to the quantum critical 
point: 



q^p{X) ^ iAr^«=,^ 

where, according to Ref. [53| : 



A^^ := A„ -)- - 2z ■ 



(113) 



(114) 



We remind that v is the correlation length critical expo- 
nent, which can in principle depend on the direction of 
the quench in the A space, and z is the dynamical criti- 
cal exponent. For a single-component, Aq is the scaling 
dimension of the operator V = d\H . Assuming that V 
is non-zero in the limit A — > from the scale invariance 
of the action we find that A„ = d + z — \ jv^ which using 
Eq. (|114p reproduces the correct scaling of the probability 
to excite the system Pox ^ lApL^/". Note that the diag- 
onal components of the geometric tensor give the fidelity 
susceptibilities Xf discussed earlier (see also Ref. [t^ for 
the review). These observations suggest a close relation 
between the scalings of Pcx with the quench amplitude. 
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the quantum geometric tensors and the fidelity suscep- 
tibility characterizing the static ground state properties 
of the system. A similar analysis shows the connection 
between the heat and the energy susceptibility xe (see 

Eq. gg])). 

The importance of the geometric tensor becomes clear 
in the contest of multi-parametric dynamics, when sev- 
eral parameters involved in the model can depend on 
time. In this case the integrals in the expressions for 
the physical quantities should be replaced by the line 
integrals in the parameter space. One can define a "re- 
sponse" in terms of one parameter by driving another 
parameter. The corresponding susceptibility of this re- 
sponse is described by the real parts of the non-diagonal 
matrix elements of the quantum geometric tensor. More- 
over the imaginary part of the geometric tensor is respon- 
sible for the Berry phase which may appear as a result 
of adiabatic evolution in this multi-parameter space. We 
are going to come back to these questions in the future. 



B. Integrable perturbations and ergodicity 

It is interesting to note that the Cauchy-Schwarz in- 
equality, although it gives simply an upper bound, it ap- 
pears to be exact for the scaling functions we derived 
here. On the other hand for some type of perturbations 
away from the QCP's (the so-called integrable perturba- 
tions) it seems possible to formulate a lower bound as 
well. It is known [s^ that for the models which have 
integrals of motion, the time averages of the correlation 
functions have a lower bound and therefore their dynam- 
ics can be non-ergodic. In particular, if the quantities 
Ai{i = 1, • • • m) are the integrals of motion, [H, Ai] = 
Vi, defined such that {Ai) = 0, then the dynamical cor- 
relation function of the operator X{t) is bounded from 
below as: 



bound for the time average of the energy susceptibility 



lim — 



{X{t)X{Q))dt > (XA) ■ (AA)-i • (XA) 



(115) 

where (AA) is the m x m matrix with elements (AiAj). 
Provided that this matrix can be diagonalized by some 
unitary transformation and that the Hamiltonian is iden- 
tified with one of the operators Ai, say H = Ai, we ob- 
tain 

the transformed operator, and therefore [& 



6ij , where the notation A stands for 



lim 1 / {xmmdt >j2^ 



> 



{XHY^ 



- {AD - m ■ 

(116) 

We implicitly assume that our Hamiltonian is defined in 
such a way that {H) = 0. 

The sine- Gordon model represents an integrable de- 
formation (perturbation) of the conformal field theory. 
Therefore it has many (infinite number) integrals of mo- 
tion. Taking a: (t) = /^^ cos(/3(/>(a;,t)) we obtain a lower 



y'^p — lim - , 

T^oo T Jo 



dx G(x, T)di 



(117) 



By the Mazur inequality (|116p this susceptibility has a 
lower bound, which is nonzero. I.e. xe can either be 
finite or diverge. The Mazur inequalities above can be 
also directly generalized for other time averages of the 
susceptibilities like xf' defined in Eq. ([55)) as well as 
higher order adiabatic susceptibilities [83 |. 



VIII. CONCLUSIONS 

In this work we analyzed the quench dynamics start- 
ing (ending) at the quantum critical point focusing on the 
sine-Gordon model. We derived the universal scalings of 
such quantities as the probability of exciting the system 
(Pex), density of excited quasiparticles (riex), entropy and 
heat densities {Sd and Q) with the quench amplitude. 
These scalings are fully determined by the critical expo- 
nents z and v characterizing the quantum critical point 
and agree with general expectation presented in [45|. 

In particular for the type of quenches where the tun- 
ing parameter changes as a power law near the quantum 
critical point: \(t) « vf^ /r\, we showed that the scalings 
of Pcxj 'T-ex, and Sd are associated with the singularities 
of generalized adiabatic susceptibilities X2r+2(A) of order 
2r-|-2 (see Eq. ([7])), while if the quench ends at the critical 
point the scaling of Q is associated with the singularity 
of X2r+i- We note that for r > it is sufficient to have 
gapless systems, not necessarily quantum critical point, 
in order to observe these singularities. 

For the two limits where the elementary excitations in 
the system are characterized by free bosons or fermions 
(for the sine-Gordon model these limits correspond to 
a particular choice of /?, see Eq. ([121)) we generalized 
our results for the finite temperature quenches. We 
showed that the structure of the singularities remains 
the same except that in the scaling relations for the den- 
sity of quasiparticles and heat one needs to substitute 
the dimensionality d ^ d — z for bosons (simultane- 
ously multiplying Q, riox by the temperature) and for 
d — > + z for fermions (simultaneously dividing Q, 7icx 
by the temperature). This changes are a direct mani- 
festation of the statistics of the quasiparticles: bunch- 
ing of bosons enhancing non-adiabatic effects and anti- 
bunching of fermions suppressing the transitions. 

We believe our results can be directly tested both nu- 
merically and experimentally, especially in setups realiz- 
able with cold atoms. 
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